Geographical effects on cascading breakdowns of scale-free networks 
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Cascading breakdowns of real networks are severe accidents in recent years, such as the blackouts 
of the power transportation networks in North America. In this paper, we study the effects of 
geographical structure on the cascading phenomena of load-carried scale-free networks, find that 
more geographically constrained networks tend to have larger cascading breakdowns. Explanations 
by the effects of circles and large betweenness of small degree nodes are provided. 
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Recently dynamical processes on networks has been 
highly concerned and widely investigated [H 0, Q 
Among many of the dynamical features of networks, ro- 
bustness attracts much attention 0, IE IE 07 much of 
which focus on scale- free (SF) networks, i.e., the degrees 
of nodes satisfy a power law distribution: P(k) ~ k~ x , 
for their ubiquity in real systems Q. The heterogeneity 
of the degrees often makes the scale-free networks sen- 
sitive to intentional attackJE 0, while it is resilience 
to random breakdowns [E 0j7 and also resilience under 
avalanche phenomena by the role of the hubs that sus- 
tain large amounts of grain, playing the role of reservoirs 
0- Furthermore, for cascading failures, the load-carried 
SF network is fragile even when one attacks only one 
node or very few nodes with the largest degrees |10| . 

Since many real networks exist in two or three dimen- 
sional physical spaces, it is helpful to study the geograph- 
ical complex networks and it has attracted much atten- 
tion recently El El El El E3 

It has been shown that 
geographical structure has great influence on percolation 
thresholds [l5|. Since many real systems bear cascading 
failures, such as power grid networks, traffic lines, In- 
ternet, etc., and also lay on the two dimensional global 
surface, the influence of geographical structures on cas- 
cading breakdowns is of highly importance and up to now 
is rarely studied. 

In this paper, we study the effects of geographical 
structure on the cascading phenomena of load-carried 
scale-free networks, in which each node carries a cer- 
tain type of load, such as power, traffic, etc., and if the 
node is broken down, its load will be redistributed to 
its neighbors. We investi gat e the Bak-Tang-Wiesenfeld 
(BTW) sandpile model HH]! 

as a prototypical model on 
a weighted lattice embedded SF (WLESF) network (ll| : 
and further study the betweenness distribution. Both 
validate that the more spatially loosely connected net- 
work is more robust under cascading failures, i.e., they 
have less huge avalanche events. The network is gen- 
erated as follows EH- It begins with an L x L lattice, 
with periodical boundary conditions, and for each node 
assigned a degree k drawn from the prescribed SF degree 
distribution P(k) ~ k~ x , k m. Then a node i is picked 
out randomly, according to a Gaussian weight function 



fi(r) = De \ a V>hJ | it selects other nodes and estab- 
lishes connections until its degree quota fcj is filled or until 
it has tried many enough times, avoiding duplicate con- 
nections. The process is carried out for all the nodes in 
the lattice. The clustering parameter A controls the spa- 
tial denseness of the connections. For large A limits, e.g. 
Ay/m 3> L, the weight function will be trivial, and the 
network becomes a SF random (SFR) network, i.e., ran- 
dom otherwise than SF degree distribution |17| . To com- 
pare, we also investigated lattice embedded SF (LESF) 
networks with nearest neighbor connections Here, 
we assume that the time scales governing the dynamics 
are much smaller than that characterizing the network 
evolvement, thus the static geographical network models 
are suitable for discussing the problems. 

The rules we adopted for sandpile dynamics are as fol- 
lows: (i) At each time step, a grain is added at a ran- 
domly chosen node i. (ii) If the height at the node i 
reaches or exceeds a prescribed threshold Z\ — hi, the 
degree of the node i, then it becomes unstable and all 
the grains at the node topple to its adjacent nodes: 
hi = hi — hf, and for each i's neighbor j: hj = hj + 1; 
during the transfer, there is a small fraction / of grains 
being lost, which plays the role of sinks without which 
the system becomes overloaded in the end. (iii) If this 
toppling causes any of the adjacent nodes to be unstable, 
subsequent topplings follow on those nodes in parallel un- 
til there is no unstable node left, forming an avalanche 
event, (iv) Repeat (i) -(iii). 

The main feature of the BTW sandpile model on Eu- 
clidean space is the emergence of a power law with 
an exponential cutoff in the avalanche size distribution, 
p(s) ~ s~ T e -s / Sc , where s is the avalanche size, i.e., the 
number of toppling nodes in an avalanche event, and s c is 
its characteristic size. In our studies, nodes toppled more 
than once in an avalanche event is seldom |9(, unless for 
the very large avalanches, which have already exceeded 
the exponential cutoffs. Thus we study the avalanche 
area, which is the number of distinct nodes that toppled 
in an avalanche event, instead of avalanche size. The 
avalanche area distribution follows the same form as that 
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FIG. 1: Number of avalanches with size j or area j, for LESF 
networks out of 10 6 avalanche events on one network config- 
uration, m = 4, N = 10 5 . 



of avalanche sizes 



p(a) 



-a/a c 
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where a is the avalanche area, and a c its characteristic 
size. A typical example is shown in Fig. ^ 

For BTW sandpile model on SFR networks, K. S. Goh 
et al. have shown that the avalanche area exponent r 
increases as A decreases, caused by the increasing num- 
ber of hubs playing the role of reservoirs. Here, we will 
demonstrate that for the densely connected scale-free ge- 
ographical networks, the reservoir effect is weakened, and 
the network has a smaller r. 

Figure [21 represents the avalanche area distribution 
for different A of LESF networks and WLESF networks 
with A = 1. It shows that as A decreases, the curve 
of avalanche area distribution is steeper, corresponds to 
larger r. These are the same as the results in Ref. [9|. 
The avalanche area exponent r for these data are fitted 
by formula ^] and is presented in Fig. |3| together with 
that of SFR networks for comparison. The data for SFR 
networks we obtained is consistent with that of 0. For 
large A large N limits, the SFR network tends to ER 
random graphs, for which r ~ 1.5 0^3; while LESF 
network tends to a super lattice, with each node has m 
neighbors; since in our studies m = 4, the network limits 
to a normal 2D lattice, which has a value of 1.01(2) for 
r, consistent with the previous results (ltllTsj. 

The avalanche area exponent for different A of WLESF 
network is shown in Fig. As A goes larger, avalanche 
area exponent r increases, the curves of avalanche area 
distribution become sharper in the double-log plot (see 
inset of Fig. 0J, which corresponds to fewer large 
avalanche events. This transition in r illuminates that 
when the network is geographically more loosely con- 
nected, it will be harder for large cascading events to 
occur. 



FIG. 2: Avalanche area distribution for LESF (left panel) and 
WLESF A — 1 (right panel) networks. For both panels, from 
up to down A = 10.0, 5.0, 4.0, 3.5, 3.0, 2.8, 2.6, 2.4. The 
loosing probability is / = 0.001, and m = 4, N = 10 5 . 10 
network realizations are carried out and for each 10 6 avalanche 
events are recorded for statistics. 



FIG. 3: Avalanche area exponent r vs the SF degree exponent 
A. The data are fitted by formula from the data presented 
in Fig. |5|and that of SFR networks. The network parameters 
and the statistics for SFR network are the same as that in 
Fig. 13 



The range of an edge is the length of the shortest paths 
between the nodes it connected in the absence of itself 
0, H3- If an edge's range is I, it will probably belong 
to an I + 1 circle. Thus the distribution of range in 
a network sketches the distribution of circles. The in- 
set of Fig. [5] shows that when the spatial constrains is 
slighter, as A goes larger, the range distribution drifts to 
larger ranges. It means that spatially loosely connected 
networks have fewer small order circles but more higher 
order circles. If there are many small order circles, the 
toppling grains are easier to meet, and the nodes with 
much less grains, i.e., fewer than z — 1, especially those 
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FIG. 4: Avalanche area exponent t vs the clusterness param- 
eter A, for A = 3.0 (squares), 5.0 (circles) and 10.0 (triangles). 
The data are fitted by formula U Inset: Avalanche area dis- 
tribution for A = 3.0, from top to bottom are LESF, WLESF 
A = 1, A = 2, and SFR networks. Dynamical and network 
parameters are the same as that in Fig. [2] 



with z — 2 or z — 3 grains, could also reach the toppling 
threshold z and topple. Larger order circles contribute 
less to this effect. The main frame of Fig. [S] shows the 
fraction of nodes toppled in avalanches that have pre- 
cisely z — 1 grains. As the network is less geographically 
constrained and has fewer small order circles, the fraction 
of toppling nodes with z — 1 grains increases, justifies our 
reasoning. This effect contributes to the large avalanche 
events of the densely connected networks, and explains 
the decrease of avalanche area exponent r as the network 
is more geographically constrained. 

In the following section, we studied the betweenncss 
distribution of these geographical networks. The be- 
tweenness, or betweenness centrality, of node i is defined 
as the total number of shortest paths between pairs of 
nodes that pass through i |2l| • If a pair of nodes has two 
shortest paths, the nodes along those paths are given 
a betweenness of 1/2 each. The betweenness distribu- 
tion for SF networks is reported to follow a power law 
Pb ~ b~ s , and for 2 < A < 3, the exponent is S w 2.2(1) 
[3 ■ We find that the betweenness distribution of LESF 
network decays much slower than that of SFR networks, 
as Fig. demonstrates for a particular case. The distri- 
butions for WLESF networks lay between them, but do 
not appear in the graph for clearness. The same holds 
for other A and m values. Thus there are more large 
betweenness nodes in LESF networks than in SFR net- 
works. To comprehend this, we plot the betweenness vs 
node's degree in Fig. For LESF networks the be- 
tweenness of the same degree is distributed much more 
diffusively, and on average are larger. It could be seen 
that even nodes with small degree k could have unusually 
large betweenness. 

When an avalanche occurs, the front of toppling nodes 



FIG. 5: Fraction of nodes that toppled after receiving only one 
grain in an avalanche event vs avalanche area. From bottom 
to top is LESF (squares), WLESF A = 1 (circles), A = 2 
(up triangles), A — 3 (down triangles), A = 5 (diamonds), 
and SFR network (left triangles). Each has 10 6 avalanche 
records on one network for statistics. A = 3, m = 4, N = 10 5 . 
Inset: Range distribution of the same networks; same symbols 
represent same networks as that in the main frame. 



FIG. 6: Betweenness distribution of the networks. A = 3.0, 
m = 2, and network size N = 10 4 , each has been averaged 
over 100 configurations. 



spread along geodesies, i.e., along the shortest paths be- 
tween nodes. Since the betweenness of a node is the num- 
ber of shortest paths passing through it, larger between- 
ness means that it will have higher possibility to receive 
grains in avalanching processes. In the above sandpile 
model, the toppling threshold is the node's degree, thus 
the node that has large betweenness but small degree 
will be easier to topple. As Fig. [7| shows, LESF net- 
work have more such nodes than SFR networks, and the 
situation changes continuously for WLESF network with 
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FIG. 7: Betweenness b vs degree k of nodes. Data are the 
same as that in Fig. [5] 

increasing A. This could also account for the decreas- 
ing avalanche area component r as the network is more 
geographically constrained. 



In conclusion, by studying avalanching processes on ge- 
ographical SF networks, we find that besides the reservoir 
effects of the hubs in SF networks, geography has great 
influences on the critical exponents of these systems. The 
decreasing avalanche area component r for the more geo- 
graphically constrained network hints high risks for such 
network to breakdown through cascading failures, since 
they have a much higher possibility to experience huge 
avalanche events, due to the denser connections and huge 
number of smaller order circles and larger betweenness of 
small degree nodes. Since many real networks that car- 
ried some kinds of loads, i.e., power, traffic, data packets, 
etc., are imbedded in the 2D global surface and highly 
clustered, our results indicate that they will suffer more 
severe risks under node failures. 
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